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Abstract. I review recent developements in lattice QCD at finite temper- 
ature, including the determination of the transition temperature T c , equation 
of state and diffenct static screening lengths. The lattice data suggest that 
at temperatures above 1.5T C the quark gluon plasma can be considered as gas 
consisting of quarks and gluons. 



Keywords: quark gluon plasma, lattice QCD 
PACS: 11.15.Ha, ll.10.Wx, 12.38.Mh, 25.75.Nq 



1. Introduction 
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One of the most challenging question in particle and nuclear physics is the one con- 
cerning the properties of strongly interacting mater at extremely high temperatures 
f— ^ and densities. We expect that at sufficiently high temperatures and densities the 

strongly interacting matter undergoes a transition to a new state, where quarks and 
gluons arc no longer confined in hadrons, and which is therefore often referred to 
as a deconfined phase or Quark Gluon Plasma (QGP). We would like to know at 



which temperature the transition takes place and what is the nature of the transition 
as well the properties of the deconfined phase, equation of state, static screening 
lengths, transport properties etc. Lattice QCD can provide first principle calcula- 
tion of the transition temperature, equation of state and static screening lengths 
(see Ref. [[TJ [2]) for recent reviews. Calculation of transport coefficients remains 
an open challenge for lattice QCD (see discussion in Refs. [[3 0]). 

One of the most interesting question for the lattice is the question about the 
nature of the finite temperature transition and the value of the temperature T c where 
it takes place. For very heavy quarks we have a 1st order deconfining transition. 
In the case of QCD with three degenerate flavors of quarks we expect a 1st order 
chiral transition for sufficiently small quark masses. In other cases there is no 
true phase transition but just a rapid crossover. Lattice simulations of 3 flavor 
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QCD with improved staggered quarks (p4) using iV r = 4 lattices indicate that the 
transition is first order only for very small quark masses, corresponding to pseudo- 
scalar meson masses of about 60 MeV [ [TU] . A recent study of the transition using 
effective models of QCD resulted in a similar estimate for the boundary in the quark 
mass plane, where the transition is 1st order [[§]. This makes it unlikely that for the 
interesting case of one heavier strange quark and two light u, d quarks, corresponding 
to 140 MeV pion, the transition is 1st order. However, calculations with unimproved 
staggered quarks suggest that the transition is 1st order for pseudo-scalar meson 
mass of about 300 MeV [ [TT] . Thus the effect of the improvement is significant and 
we may expect that the improvement of flavor symmetry, which is broken in the 
staggered formulation, is very important. But even when using improved staggered 
fermions it is necessary to do the calculations at several lattice spacings in order 
to establish the continuum limit. Recently, extensive calculations have been done 
to clarify the nature of the transition in the 2+1 flavor QCD for physical quark 
masses using N t — 4, 6, 8 and 10 lattices. These calculations were done using the 
so-called stout improved staggered fermion formulations which is even superior to 
other more commonly used improved staggered actions (p4, asqtad) in terms of 
improvement of flavor symmetry. The result of this study was that the transition 
is not a true phase transition but only a rapid crossover [ I12j . Even-though there 
is no true phase transition in QCD thermodynamic observables change rapidly in 
a small temperature interval and the value of the transition temperature plays an 
important role. The flavor and quark mass dependence of many thermodynamic 
quantities is largely determined by the flavor and quark mass dependence of T c . 
For example, the pressure normalized by its ideal gas value for pure gauge theory, 
2 flavor, 2+1 flavor and 3 flavor QCD shows almost universal behavior as function 
ofT/T c [|9]. 

The chiral condensate (ipip) and the expectation value of the Polyakov loop (L) 
are order parameters in the limit of vanishing and infinite quark masses respectively. 
However, also for finite values of the quark masses they show a rapid change in 
vicinity of the transition temperature. In Figure [T] I show the chiral condensate 
and the Polaykov loop as function of the temperature calculated for the pA action 
and light quark mass mi = 0.1m s , with m s being the physical strange quark mass. 
Note that in Figure we show the renormalized Polyakov loop defined as in Ref. [ 
[13] . The details of the calculations calculated can be found in Ref. [ [14] . We see 
that the chiral condensate and the renormalized Polyakov loop show rapid change 
at T c suggesting that the chiral and the deconfincmcnt transitions happen at the 
same temperature. To determine the value of the transition temperature and to 
study the interplay between the chiral and the deconfinement transition one usually 
calculates the disconnected part of the chiral susceptibility and the Polyakov loop 
susceptibility defined as 

^± = Num?) - <v^) 2 ), H = mm - (i) 2 ) (i) 



as function of the of the bare gauge coupling (3 = Q/g 2 - N a is the spatial size 
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Fig. 1. The renormalized Polyakov loop L ren (T) [[14] and the chiral condensate 
normalized to the zero temperature chiral condensate [ [TS] as function of the tem- 
perature calculated on 16 3 x 4 lattices. 



of the lattice. The susceptibilities have a peak at some pseudo-critical coupling 
f3 c . The chiral and the Polyakov loop susceptibility have been studied using lattice 
with temporal extent N T = 4 and N T = 6 and several values of the light quark 
masses m; = 0.05to s , 0.1m s , Q.2m s and 0Am s [[7J. Note that the smallest value 
of mi correspond to pion masses of about 140MeV. We find that within accuracy 
of the calculations pseudo-critical couplings (3 C determined from the disconnected 
part of the chiral susceptibility and the Polyakov loop susceptibility coincide. This 
again shows that the chiral and the deconfinement transition happen at the same 
temperature. To determine the transition temperature we have to calculate the 
lattice spacing in terms of some physical quantity. In the past the string tension 
has been used to set the lattice spacing. A more accurate determination of the 
lattice spacing is provided by the so-called Sommer scale ro defined from the static 
quark anti-quark potential as r 2 \ r =r = 1-65. Analysis of the quarkonium 
spectroscopy on the lattice lead to the value r = 0.469(7)fm [ 16j. In figure [5] I 
show the transition temperature in units of ro for different quark masses [[7J and two 
different lattice spacings, corresponding to N T = 4 and N T = 6 lattices. Note that 
the value of T c calculated at two different lattice spacings are clearly different. The 
thin error-bars in Figure [2] represent the error in the determination of the lattice 
spacing a, i.e. the error in ro/a. There is also an error in the determination of 
the gauge coupling constant C — 6/g 2 . The combined error is shown in Fig. lb 
as a thick error-bar. For Nt — 4 calculations the error is dominated by the error 
in lattice spacing, while for Nt = 6 it is dominated by the error in f3 c . With the 
data on tqT c a chiral and continuum extrapolation has been attempted using the 
most simple Ansatz r T c {ra Vl N t ) = r T c \^ t al + A(r m 7T ) d + B/N 2 . From this 
extrapolation on gets the continuum value T c r = 0.457(7) [+8] [—2] for the physical 
pion mass m^rQ = 0.321 [ [7]. The central value was obtained using d = 1.08 
expected from 0(4) scaling. To test the sensitivity to the chiral extrapolations 
d = 2 and 1 have also been used. The resulting uncertainty is shown as second 
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Fig. 2. The transition temperature in units of the tqT c from Ref. [[7] as function 
of the pion mass. 

and third error in square brackets. Using the best know value of ro — 0.469(7)fm 
we obtain T c = 192(7) (4) MeV which is higher than the most of the previous 
values. It is also significantly higher than the chemical freezout temperature at 
RHIC [ Q7] . Note that the large value of the transition temperature is mostly due 
to the large value of the string tension which is related to the Sommer scale as 
r o\/v = 1-114(4) [US]. Using the above value of r we get ^fa — 468 MeV which 
is more than 10% larger than the value yfa = 420 MeV which was used in Ref. 
[ [19] and let to T c ~ 173(8)MeV. Recently the transition temperature has been 
determined using the so-called stout staggered action and N t = 4, 6, 8 and 10 [ 
[20] . The deconfinement temperature has been found to be 176(3) (4) MeV [ 120] . 
The central value is considerably smaller than the one obtained with pA action but 
taking into account the errors the deviation is not very significant. Calculations on 
Nt — 8 lattices with pA action are needed to clarify this issue. The authors of Ref. [ 
[20] use a different definition of the chiral susceptibility which resulted in the chiral 
transition temperature of T c hi ra i — 151(3)(3)MeV. Using the chiral susceptibility 
defined above would result in a larger value of the transition temperature. 

Lattice calculations of equation of state were started some twenty years ago. 
In the case of QCD without dynamical quarks the problem has been solved, i.e. 
the equation of state has been calculated in the continuum limit [ 21j. At tem- 
peratures of about 4T C the deviation from the ideal gas value is only about 15% 
suggesting that quark gluon plasma at this temperate is weakly interacting. Per- 
turbative expansion of the pressure, however, showed very poor convergence at this 
temperature [H2]. Only through the use of new resummed perturbative techniques 
it was possible to get agreement with the lattice data [[23] [24] [25]. To get a reliable 
calculation of the pressure and the energy density improved action have to be used 
[ [26] [27] . Very recently calculations with the so-called asqtad and pA action have 
been done on lattices with temporal effects N T = 4 and 6 [ [28], [29] . In Figure [3] 
the interaction measure e — ip is shown as function of the temperature for the pA 
action. Calculations performed for N T — 4 and N T — 6 give similar results. This 
means that cutoff effects arc under control. Furthermore, there is a good agreement 
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Fig. 3. The interaction measure calculated for the p4 action [[29]. 

between pA and aqstad calculations. We see that close to T c the interaction measure 
is very large, which means that quark gluon plasma at this temperature is very far 
from the conformal limit. At high temperature the value of the interaction measure 
is consistent with the perturbative estimate. 

2. Spatial correlation functions 

To get further insight into properties of the quark gluon plasma one can study 
different spatial correlation functions. One of the most prominent feature of the 
quark gluon plasma is the presence of chromoelectric (Debye) screening. The easiest 
way to study chromoelectric screening is to calculate the singlet free energy of static 
quark anti-quark pair (for recent review on this see Ref. [I30j). which is expressed 
in term of correlation function of temporal Wilson lines 

cxp(-Fi (r, T)/T) = Tv{W {r)W^ (0)) . (2) 

L = TrW is the Polyakov loop. In absence of dynamical quarks the free energy grows 
linearly with the separation between the heavy quark and anti-quark in the confined 
phase. In presence of dynamical quarks the free energy is saturated at some finite 
value at distances of about 1 fm due to string breaking [[30] [3T]. Above the decon- 
finement temperature the singlet free energy is exponentially screened, at sufficiently 
large distances [[3J, i.e. i*\ (r, T) = F^T) - exp(-m D (T)r). The inverse 

screening length or equivalently the Debye screening mass mo is proportional to the 
temperature. In leading order of perturbation theory it is mj = y/l + Nf/3g(T)T. 
Beyond leading order it is sensitive to the non-perturbative dynamics of the static 
chromomagnetic fields. The Debye screening mass has been calculated in pure 
gauge theory (Nf = 0) [ [32] and in 2 flavor QCD (Nf — 2) [ [33] and is shown 
in Figure^] at different temperatures. The temperature dependence of the lattice 
data have been fitted with the simple Ansatz motivated by the leading order re- 
sult : mu(T) = Ay/l + Nf /3g(T)T. Here g(T) is the two loop running coupling 
constant. This simple form can fit the data very well and we get A ~ 1.4 both 
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Fig. 4. The Debye mass calculated in quenched QCD [ [32] and 2 flavor QCD [ 
I3"5] . The lines show the leading order fit together with its uncertainty. 

for Nf — and Nf = 2. Thus the temperature dependence as well as the flavor 
dependence of the Debye mass is given by perturbation theory. We also see that 
non-perturbative effects due to static magnetic fields significantly effect the electric 
screening, resulting in about 40% corrections. However, the non-perturbative cor- 
rection is the same in full QCD and pure gauge theory. Let us note that in SU(2) 
gluodynamics the corrections to the Debye mass are even larger, the Debye mass 
is 1.6 times larger than the leading order result [(M[ [3H1 EH]- This situation can 
be understood in terms of dimensionally reduced effective theory, where the effect 
of hard modes with momentum p ~ irT is integrated out and which contain only 
static electric and magnetic fields [ [37] . The validity of dimensional reduction has 
been tested in a wide temperature range [ (3S] [3S] . 

At zero temperature the static quark anti-quark potential is determined from 
the Wilson loops: V(r) — — 1/t In W(r, t), t — > oo. At large separation the Wilson 
loop obeys the area law W(r, t) ~ exp(— art) which means that the potential grows 
linearly with distance r. At finite temperature we can consider the spatial Wilson 
loops. They obey area law at any temperature W s (x, z) ~ exp(— a s (T)xz) [[38] [39], 
Below the transition temperature the spatial string tension is very close to the 
usual zero temperature string tension. Well above the deconfinement transition 
temperature the spatial string tension is expected to be y/a s (T) = cm9 2 {T)T [[39]. 
This is because in dimensionally reduced theory it is given by Cm<?§ and at leading 
order the 3-dimensional gauge coupling is g\ = g 2 (T)T. The spatial string tension 
has been calculated on the lattice in quenched QCD (Nf = 0) [[3D] and 2+1 flavor 
QCD [ [¥T] and the results are shown in Figure [5] The lattice data can be fitted 
very well with the simple form : y/a s (T) — cm9 2 {T)T. Here again g(T) is the 
2-loop running coupling. For the fit we get the value of cm which agrees well with 
the result of dimensional reduction [ [41] . The 3-dimensional gauge coupling g\ has 
been calculated more systematically in perturbation theory and also led to a very 
good agreement with the lattice data [ [35] . 
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Fig. 5. The spatial string tension calculated in quenched QCD [[40] and 2 flavor 
QCD [ 41 . The lines show the leading order fits. 

3. Conclusions 

In recent years significant progress has been made in calculating bulk thermody- 
namic observables on the lattice as well as spatial correlation functions. This calcu- 
lations suggest that at temperatures T > 1.5T C thermodynamics can be described 
reasonably well using weak coupling approaches: resummed perturbation theory and 
dimensional reduction. The temperature and flavor dependence of static screening 
length is well described by perturbation theory. However, the value of the screen- 
ing lengths is to large extent non-perturbative and influenced or determined by the 
dynamics of static magnetic fields. Furthermore, there is no evidence for the large 
value of the gauge coupling constant at scale T. Clearly more precise lattice data 
and further perturbative calculations are needed to establish the nature of quark 
gluon plasma in the temperature, T > 1.5T C . 
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